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Abstract 

We consider proper 5-colourings of the kagome lattice. Proper (/-colourings cor- 
respond to configurations in the zero-temperature (/-state anti-ferromagnetic Potts 
model. Salas and Sokal have given a computer assisted proof of strong spatial 
mixing on the kagome lattice for q > 6 under any temperature, including zero tem- 
perature. It is believed that there is strong spatial mixing for q > 4. Here we give 
a computer assisted proof of strong spatial mixing for q = 5 and zero temperature. 
It is commonly known that strong spatial mixing implies that there is a unique 
infinite-volume Gibbs measure and that the Glauber dynamics is rapidly mixing. 
We give a proof of rapid mixing of the Glauber dynamics on any finite subset of the 
vertices of the kagome lattice, provided that the boundary is free (not coloured). 
The Glauber dynamics is not necessarily irreducible if the boundary is chosen arbi- 
trarily for q = 5 colours. The Glauber dynamics can be used to uniformly sample 
proper 5-colourings. Thus, a consequence of rapidly mixing Glauber dynamics is 
that there is fully polynomial randomised approximation scheme for counting the 
number of proper 5-colourings. 

1 Introduction 

Proper colourings correspond to configurations in the zero-temperature anti-ferromagnetic 
Potts model. In this paper we will show that the system specified by proper 5-colourings 
of the kagome lattice has strong spatial mixing, and that the Glauber dynamics is rapidly 
mixing. The previously best known result [20] on mixing on the kagome lattice was for 6 
colours. It is believed [20] that there is strong spatial mixing for 4 or more colours, and 
hence our result is narrowing the gap between what is believed and known. In Section 1.1 
below we give an introduction to mixing, and in Section 2 we state our results and discuss 
related work. 

*Some of the work has been done at the Department of Computer Science, University of Warwick, 
Coventry, CV4 7AL, UK. 
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Figure 1: (a) The kagome lattice, here drawn in a coordinate system. We illustrate graphs 
such that a face represents a vertex, (b) A vertex (x, y) G V Q dd and its right neighbour, 
(c) A vertex (x, y) G V eYeQ and its four neighbours. 



1.1 Definitions and background 

The kagome lattice, Figure 1(a), is a natural lattice of interest in statistical physics [20]. 
Instead of drawing graphs in the traditional way, with a vertex denoted with a solid 
circle and an edge denoted with a line segment, we draw graphs such that faces represent 
vertices. Two adjacent faces therefore represent two adjacent vertices. Let Q = (Vg : Eg) 
denote the kagome lattice with vertex set Vg and edge set Eg. We have Vg = V oAd U V cven , 
where 

V dd = {( x ,y) | x, y G Z are both odd}, 

Kven = {{x, y) I x = 4fci +r, y = Ak 2 + r, k u k 2 G Z, r G {0,2}}. 

The edge set 

Eg = {((x,y),(x + 2,y)) | (x,y)eV odd } U 

{((x,y),(x + k 1 ,y + k 2 )) | (x, y) G V cvcn , h G {-1, 1}, k 2 G {-1, 1}}. 

Note that both vertices in Figure 1(b) are in V a dd, and we see that two adjacent vertices 
in V dd differ by 2 in their x-coordinate. In Figure 1(c), the centre vertex is in V even , and 
we see that its four neighbours are in V dd- 

A region R C Vg is a finite non-empty subset of the vertex set of the kagome lattice. 
The subset OR C Vg denotes the vertex boundary of R such that OR is the set of vertices 
that are not in R but are adjacent to any vertex in R. The edge set E(R) is the set of all 
edges (u, v) G Eg such that at least one of the vertices u and v is in R. The edge boundary 
£R of R is the set of all edges (w, v) G E(R) such that exactly one of the vertices u and 
v is in it! and the other one is in OR. 

The set Q — {1, . . . , q} denotes the set of q colours, and the set Q$ = {0} U Q. The 
colour represents "no colour". A q- colouring of a region R is a function from R to the 
set Q, and a go-colouring of R is a function from R to Qo- A 0-colouring of R is a function 
from R to the set {0}, which means that all vertices in R are assigned colour 0. We often 
write only colouring when it is obvious from the context if it is a q-, q$- or 0-colouring, 
or if any colouring will do. Let o be a colouring of a region R. If R' is a subset of R 
then a(R') is the colouring of R' induced by a. Furthermore, for a vertex v G R, a{y) 
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is the colour of v under a. Let fl~^ denote the set of all g-colourings of the region R. 
For two colourings a, a' G the Hamming distance between a and a' is the number 
of vertices in R on which a and a' differ. A colouring a of R is proper if no adjacent 
vertices receive the same colour. That is, a{u) ^ a(v) for all adjacent vertices u and v in 
R. Let Qr denote the set of all proper g-colourings of the region R. Given a go-colouring 
B of dR, a proper g-colouring a of R agrees with B if a(u) ^ B(v) for all (u, v) G ER, 
where u G R. We let Qr(B) denote the set of all proper g-colourings of R that agree with 
B. The uniform distribution on Qr(B) is denoted ir^, and for any subregion R' C R, let 
t^b,r' denote the distribution on proper (/-colourings of R' induced by n B . 

In this paper we will show that the system specified by proper 5-colourings of the 
kagome lattice has strong spatial mixing. Informally, strong spatial mixing means that if 
R is a region and B is a go-colouring of dR, then the effect the colour of a vertex w G dR 
has on a vertex v G R decays exponentially with the distance between w and v. The 
effect is measured with the total variation distance. For two distribution D\ and D 2 on 
a set S, the total variation distance between D\ and D 2 is defined as 

The following definition of strong spatial mixing is taken from [12] and is adapted to the 
kagome lattice. 

Definition 1 (Strong spatial mixing). The system specified by proper g-colourings of the 
kagome lattice has strong spatial mixing if there are two constants a > and e G (0, 1) 
such that, for any region R, any subregion R' C R, any two go-colourings B and B' of dR 
which differ on exactly one vertex w G dR and such that B(w) ^ and B'(w) ^ 0, 

dTv(vr B ^,vr B , R < a\R'\(l - e) d ^ R '\ 

where d(w, R') is the minimal distance within R from w to some vertex of R' . 

A distribution 7r on the set of proper g-colourings of the infinite kagome lattice is 
an infinite-volume Gibbs distribution if, for any region R and any proper g-co louring a 
of the kagome lattice, the conditional distribution ir(-\a(Vg\R)) on flu (conditioned on 
the colouring a(Vg\R) of all vertices other than those in R) is 7r B , where B = a(dR). 
It is known that there is always at least one infinite-volume Gibbs distribution, and the 
question of interest is to determine whether it is unique or not. This question is central 
in statistical physics because it corresponds to the number of macroscopic equilibria for a 
given system. The phenomenon of non-uniqueness corresponds to what is referred to as a 
phase transition. A consequence of strong spatial mixing is that the infinite- volume Gibbs 
distribution is unique [8, 23, 24]. For more on Gibbs distributions, see for example [9] 
or [10]. 

Another question of interest is to determine how quickly the system converges to 
equilibrium. The answer to this question is connected to the quantities a and e in 
Definition 1 above. From a statistical physics point of view, this question is important 
for understanding phenomena such as how the system returns to equilibrium after a 
shock forces it out of it. In this paper we consider a famous dynamical process called 
the Glauber dynamics which models how the system converges. The Glauber dynamics, 
defined next, is a Markov chain that performs single-vertex heat-bath updates. 
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Figure 2: A 2- vertex region of the kagome lattice with different colourings. A vertex is 
labelled with its colour. The two colourings in (a) are "frozen" and do not communicate 
in the Glauber dynamics for q — 5. However, when restricting the boundary to the 
0-colouring in (b), the two colourings do communicate. 

Definition 2 (Glauber dynamics). For any region R and any (? -colouring B of dR, the 
Glauber dynamics is a Markov chain with state space Qr(B), and a transition is made 
from a state a to a' in the following way: 

1. Choose a vertex v uniformly at random from R. 

2. Let Q v be the set of colours which are assigned to the neighbours of v (either in a 
or B). 

3. Choose a colour c uniformly at random from Q\Q V and obtain the new colouring 
o~' from a by assigning colour c to vertex v. 

A sufficient condition for the Glauber dynamics to be connected (that is, any proper 
colouring can be obtain from another proper colouring by a series of transitions) is to 
have q > 6. In general, with the Glauber dynamics defined similarly on any underlying 
infinite graph of maximum degree A, having q > A + 2 is a sufficient condition for 
the dynamics to be connected. In this paper we focus on 5-colourings and in order to 
guarantee that the Glauber dynamics is connected we will have to restrict the colourings 
B of the boundary to the 0-colouring (see Figure 2). For this reason, our 5-colour mixing 
result for the Glauber dynamics is restricted to the 0-colouring of the boundary. It is 
worth pointing out that if we add moves to the Glauber dynamics that allow swapping 
the colours of two neighbouring vertices (when this move is allowed with respect to the 
colouring of the rest of the vertices) then the new dynamics is connected for any region 
R and any g - c °l° urm g of dR if q > 5. This fact is true for any graph of maximum 
degree A and q > A + 1 colours. This augmented Glauber dynamics can be simulated by 
the heat-bath dynamics on edges which we define as follows: Choose an edge e = (vi,v 2 ) 
uniformly at random and simultaneously recolour v\ and i>2 uniformly at random from 
the allowed colourings. 

If the Glauber dynamics is connected, and hence ergodic, then is the unique sta- 
tionary distribution. This follows from the fact that the Glauber dynamics is reversible 
with respect to 7r B . For the same reason, njg is the unique stationary distribution of the 
heat-bath dynamics on edges. The Glauber dynamics can be used as a sampler to sample 
colourings from the uniform distribution on Qr(B). This can be done efficiently if the 
Glauber dynamics is rapidly mixing (see definition below), which means that it quickly 
reaches its stationary distribution. 

Definition 3 (Mixing time). Consider the Glauber dynamics on a region R with bound- 
ary colouring B. Let P t (a, a 1 ) be the probability of going from state a to a' in exactly t 



4 



steps. For any 5 > 0, the mixing time 

t(5) = max min{t : d'ry(P t (o-, O^b) < ^ f° r an t > t }. 

aen R (B) 

The Glauber dynamics is rapidly mixing if t(S) is upper-bounded by a polynomial in the 
region-size \R\ and log (1/5). 

It is a well-known fact that if the system has strong spatial mixing then the Glauber 
dynamics is (often) rapidly mixing [8, 17, 23]. In Section 7 we will study this fact and 
see how strong spatial mixing and rapid mixing are closely related. For q > 6 colours (or 
q > A + 2 in general) it is straightforward to apply Theorem 8 in [12] in order to infer 
rapid mixing from strong spatial mixing. However, with q = 5 colours we cannot rely 
entirely on previous results. We will establish certain properties of 5-colourings of the 
kagome lattice and show that the Glauber dynamics is rapidly mixing for q = 5 under 
the 0-colouring of the boundary. 

In [15] it is explained how approximate counting and almost uniform sampling are 
related. If there is a method for sampling (almost) uniformly at random in polynomial 
time from the set of proper colouring of a finite region R, then we can construct a fully 
polynomial randomised approximation scheme, or FPRAS, for counting the number of 
proper colourings of R. Thus, if the Glauber dynamics is rapidly mixing then we could 
use it to construct (in a non-trivial way) an FPRAS for estimating \Qr\. For details on 
the topic of how sampling and counting are related, see Jerrum [15] and Jerrum, Valiant 
and Vazirani [16]. 

2 The results and related work 

We will prove the following theorems, which improve previously known results on mixing 
for proper colourings of the kagome lattice. 

Theorem 4. The system specified by proper 5-colourings of the kagome lattice has strong 
spatial mixing. 

Theorem 5. For any region R of the kagome lattice and q = 5 colours, the Glauber 
dynamics is rapidly mixing on R under the 0-colouring of dR. The mixing time t(S) G 
0(n 2 + n log |) ; where n is the number of vertices in R. 

Theorem 6. For any region R of the kagome lattice and q = 5 colours, the heat-bath 
dynamics on edges is rapidly mixing on R under any q Q -colouring of dR. The mixing 
time t(8) G 0{n 2 + n log |) ; where n is the number of vertices in R. 

The previously best known result on mixing on the kagome lattice is that of Salas and 
Sokal [20]. They provided a computer assisted proof of strong spatial mixing for q = 6 
colours. It is believed [20] that there is strong spatial mixing for q > 4 colours. 

It is worth mentioning some previous general results on mixing. Independently, Jer- 
rum [14] and Salas and Sokal [20] proved that for proper g-colourings on a graph of 
maximum degree A the Glauber dynamics has 0(n log n)-mixing for q > 2 A, where n is 
the number of vertices of the region. For q = 2 A, Bubley and Dyer [3] showed that it 



5 



mixes in 0(n 3 ) time and Molloy [18] showed that it mixes O(nlogn) time. Vigoda [22] 
used a Markov chain that differs from the Glauber dynamics and showed that it has 
0(n logn)-mixing for q > (11/6) A. This result implied that the Glauber dynamics is 
rapidly mixing for q > (11/6)A. Goldberg, Martin and Paterson [12] showed that any 
triangle- free graph has strong spatial mixing provided q > a A — 7, where a is the solution 
to a a = e (a « 1.76322) and 7 = 4a3 ~ ( 6 ° 2 2 ~ L 3 ) Q+4 ~ 0.47031. Note that their result cannot 
be applied to the kagome lattice since its edge set contains triangles. However, for other 
4- regular graphs, such as the square lattice Z 2 , it follows that mixing occurs for q > 7 
colours. The technique Goldberg, Martin and Paterson used in [12] is well suited to be 
extended to involve special cases that depend on the particular graph under considera- 
tion. Involving such special cases can improve the mixing bounds. In order to deal with 
all special cases it might be helpful to incorporate computer assistance. This has been 
done in [12] for the lattice Z 3 . The general result gives mixing for q > 11 colours but 
by taking advantage of the geometry of the lattice it has been shown that mixing occurs 
for q > 10. This proof is computer assisted. Another computer assisted proof of mixing 
in [12] is given for the triangular lattice and q — 10 colours. This result was improved by 
Jalsenius [13] to q = 9 by exploiting the geometry of the lattice even further. Goldberg, 
Jalsenius, Martin and Paterson used the technique from [12] and gave in [11] a computer 
assisted proof of mixing for q = 6 on the square lattice Z 2 . This is an alternative proof 
of the result of Achlioptas, Molloy, Moore and van Bussel [1] (who also used computer 
assistance). In this paper we will refine the technique Goldberg, Martin and Paterson 
introduced in [12] to show mixing on the kagome lattice for q = 5 colours. Both the 
square lattice and the kagome lattice are 4-regular graphs, but the kagome lattice con- 
tains triangles whereas the square lattice does not. An interesting observation is that the 
presence of triangles seem to have a positive effect on the technique we use to show strong 
spatial mixing. Attempts to prove mixing with 5 colours on the square lattice with this 
technique has failed so far. The absence of triangles seem to be one strong reason why 
(assuming the square lattice does have strong spatial mixing with 5 colours). 

3 The framework 

When Goldberg, Martin and Paterson [12] derived improved mixing bounds for spin 
systems consisting of proper colourings, they introduced the notion of a vertex-boundary 
pair. A vertex-boundary pair is a data structure holding information about a region R 
and colourings of dR. The idea is to derive certain properties of the vertex-boundary 
pairs which can be easily translated into properties such as whether there is strong spatial 
mixing or not. When Goldberg, Martin and Paterson derived these properties, it turned 
out to be convenient to work with edge-boundary pairs. An edge-boundary pair (defined 
in the next section) contains colourings of the edge boundary ER rather than the vertex 
boundary dR. 

Definition 7 (Vertex-boundary pair). A vertex-boundary pair X consists of 

• a region R X) 

• a distinguished boundary vertex w% £ dRx, and 
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• a pair (Bx, B'x) of ^-colourings of dRx that are identical on all vertices except on 
wx, where they differ. The colour of wx is in Q for both Bx and B' x . 

Note that the colour of the distinguished vertex wx has to be in the set Q. That is, 
Bx(wx) 7^ and B' x {wx) 7^ 0. Definition 1 of strong spatial mixing can be rephrased 
using the definition of a vertex-boundary pair. That is, in order to show strong spatial 
mixing, we will show that there are two constants a > and e G (0, 1) such that for 
every vertex-boundary pair X and every subregion R! C R x , 

dTv(n Bx ,R',n B/xiR/ ) < a\R'\(l - e) d ^ R '\ 

One approach to show exponential decay of the total variation distance in the distance 
between wx and Rl is to construct a suitable coupling (defined next) of the distributions 
hb x and 7Tg^, . For two distributions D\ and D 2 on a set S, a coupling \1/ of D\ and D 2 is 
a joint distribution on S x S with marginal distributions D\ and D 2 . If the pair (X±, X 2 ) 
is a random variable drawn from \1/ then 

d TV (D 1 ,D 2 )<Pi[X 1 ^X 2 ]. 

Thus, in order to upper-bound the total variation distance, one can find some suitable 
coupling ^ and compute the probability of having X\ 7^ X 2 . The aim here is to construct 
a coupling ^x of ttb x and itb> x such that if the pair (a, a') of colourings is drawn from 
^x then the probability that er and a' differ on R' C R x decreases exponentially with 
the distance between the discrepancy vertex wx £ dRx and R'. For a vertex v G Rx we 
define the indicator random variable l-qi XtV for the event that the colour of v differs in a 
pair of colourings drawn from ^x- Hence, the quantity Yl v eR x * s the expected 

number of vertices in Rx on which the colours differ in a pair of colourings drawn from 
^ x- If * s small enough for all vertex-boundary pairs X and vertices v G Rx then 

we can infer strong spatial mixing (Section 6) and rapid mixing (Section 7). 

4 Edge discrepancies 

Similarly to the definition of a vertex colouring we define a q-, qo- and 0-colouring of a set 
E C Eg of edges to be a function from E to Q, Qo and {0}, respectively. If B is an edge 
colouring of E, and E' is a subset of E then B(E') is the colouring of £" induced by B. 
For an edge e G -E 1 , -B(e) is the colour of e under B. Given a region i? and a go- c °l° ur i n g 
.B of £R, a proper g-colouring a of R agrees with B if a(u) 7^ -B(e) for all edges e G £^-R, 
where u G -R is incident to e. We let Vt R (B) denote the set of all proper g-colourings of 
R that agree with B. The uniform distribution on VLr(B) is denoted -kb- 

Let E be a set that contains the four edges that are incident to some vertex v G Vg. 
Two edges e, e' G E are adjacent if there is a clockwise ordering around v of the edges 
in E such that e' follows immediately after e. Similarly to a vertex-boundary pair we 
define an edge-boundary pair X as follows. Note that this definition is equivalent to the 
notion of a relevant boundary-pair in [12]. 

Definition 8 (Edge-boundary pair). An edge-boundary pair X consists of 

• a region Rx, 
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• a distinguished boundary edge ex = (wx,vx) G £Rx with wx G dRx, vx G -Rx, 
and 

• a pair (B X ,B' X ) of g -colourings of £R X that are identical on all edges except on 
ex, where they differ. 

We require 

• B x (e x ) G Q and B' x (e x ) G Q, and 

• any two adjacent boundary edges that share a vertex in <9i?x have the same colour 
in at least one of the two colourings Bx and B' x (and so in both of Bx and B' x 
except when edge ex is involved). 

Suppose X is an edge-boundary pair. For a coupling ^/x of ttb x and 7Tb^ we define 
t° be the indicator random variable for the event that, when a pair of colourings 
is drawn from ^x, the colour of vertex v G Rx differs in these two colourings. For any 
edge-boundary pair X we define \l/™ in to be some coupling of ti Bx and n B i minimising 
E[1# X) „ X ]. For every pair of colours c, d G Q, let p x m (c,c') be the probability that 
C(vx) = c and C'(v x ) = c', where (C, C") is a pair of colourings drawn from \l/™ m . For 
a vertex v G let d(ex,v) denote the distance within Rx from edge ex to v. Thus, 
d{ex,vx) = 1 and if v G i?x adjoins vx then d(ex,f) = 2, and so on. We wish to 
construct a coupling ty x of 7Tb x and 7r B ^ such that E[l^ Xj „] decreases exponentially in 
the distance rf(ex, f)- In order to do this we use a recursive coupling. To aid the analysis 
we define a labelled tree T x associated with each edge-boundary pair X. The notion of 
T x was introduced by Goldberg, Martin and Paterson in [12]. 

Suppose X is an edge-boundary pair. We will now construct the tree Tx- Start with 
a node r which will be the root of T X - For every pair c, c' G Q of distinct colours, add 
an edge labelled (p™ in (c, c'),vx) from r to a new node r CiC /. Let ei, e^ be the clockwise 
ordering of the edges incident to vx (excluding edge ex) such that ex appears between 
e3 and e±. The i-th neighbour of vx denotes the vertex that is incident to e^. If the i-th 
neighbour of Vx is not in Rx then we define Xj(c, c') = 0. If the i-th neighbour of vx is 
in R x then let X^c, c') be the edge-boundary pair consisting of 

• The region R Xi (c,c>) = Rx\{v x }, 

• the distinguished boundary edge ex 4 ( c , c ') = e^, and 

• the pair (B Xi (c, c '), B'xAc &)) °^ <?o-colourings of ERx^d) such that both colourings are 
identical to B x on all edges in £R Xi (c,c')\{ei, e 2, e z\- The colours of the boundary 
edges in {e\ 1 e2 1 e^} are assigned as follows. 

— B Xl (c,c')(ei) = d and B' x ^ cl) ( ei ) = c. 

— For the boundary edge ej G {ei,e 2 ,e 3 } such that j < i, both B Xi ( c ,c')( e j) and 

— For the boundary edge G {ei,e2,es} such that j > i, both -Bxi(c,c')( e j) and 

5 X i (c,C)( e i) ale C - 
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If the i-th neighbour of Vx is in R x , recursively construct the tree T Xl (c,c') and join it to 
Tx by adding an edge with label (1, •) from r cc i to the root of TxA c ,cf)- Note that if vx 
has no neighbours in Rx then r c c i is a leaf. That completes the construction of Tx- 

We say that an edge e of Tx is degenerate if the second component of its label is "•" . 
For edges e and e' of T x , we write e — > e' to denote the fact that e is and ancestor of e'. 
That is, either e = e', or e is a proper ancestor of e'. Define the level of an edge e of Tx to 
be the number of non-degenerate edges on the path from the root down to, and including, 
e. Suppose that e is an edge of Tx with label (p, v). We say that the weight w(e) of edge 
e is p. Also the name n(e) of edge e is v. The likelihood 1(e) of e is Yle'-e'->e w ( e )- The 
cos/ 7(f,Tx) of a vertex v G i?x is E e -n(e)=i> K e )- ^ the region i?x is not connected and 
vertex vx and a vertex t> G i?x belong to different connected components, then there will 
be no edge with name v in Tx and we define 7(1;, Tx) = 0. We have the following lemma, 
which is proved in [12] as Lemma 12. 

Lemma 9 ([12, Lemma 12]). For every edge-boundary pair X there exists a coupling ^x 
of n Bx and n B i x such that E[1^ X) „] < 7(1;, Tx) for all v G Rx- 

A key ingredient from the construction of Tx that affects j(v,Tx) is the quantity 
Efl^mm ^], which we denote v(X). Thus, 

v(X)=E[l v $n tVx ]= Px n (c,c>). 

c,c'eQ, 

c^c' 

For an edge-boundary pair X and an integer d > 1, let E d (X) denote the set of 

level-d edges in T x , and define T d (X) = E eG E d (x) Z ( e )- We define = for d > 1. 

Equivalently, we can define T d (X) recursively: 

T\X) = u(X)= J2 P7 n (c,c'), (1) 
c,c'eQ, 

and for d > 1 we have 

T d (X)= J2 PT^c'^T^iX^c')). (2) 

c,c'eQ, i=l 

Lemma 10. Suppose X is an edge-boundary pair and R! C T/ien there is a coupling 
ofn Bx and 7r B t x such that 

fGR' d>d(e x ,R') 

Proof. By Lemma 9 there is a coupling ^x of 7Tb x and ttb' such that E[1^ X) „] < 7(1;, T x ) 
for t> G i?x- Thus, 

£ E[l* x ,„] < £ 7 (^x) < £ £ /(e) 

nei?' v£R' veR' e:n(e)=v 

< £ £ '(«)< £ r-(A-). 

d>d(e x ,R') eeE d (X) d>d(e x ,R') 

□ 
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5 Exponential decay of F d (X) 



Suppose X is an edge-boundary pair. Let B be the colouring of SRx such that B(e) = 
Bx(e) for e G SRx \ { e x} and B(ex) = 0. For i e Q, we define nj(X) to be the 
number of proper g-colourings a in such that a{vx) = i- For e Q, we define 

= Ejeo\{i,i'}«iW and 



n,(X) + iV M ,(X)' 



Suppose c = B x (ex) and c' = B' x (e x ). Then /i CjC /(X) is the probability that vx receives 
colour c in ~Kb' x i an d /i c ' jC (X) is the probability that vx receives colour d in tt Bx . We now 
define fx(X) = max[ / u CjC /(X), /v, c PO]- 

Lemma 11. For every edge-boundary pair X , v(X) < fi(X). 

Proof. Let X be an edge-boundary pair and suppose without loss of generality that 
Bx(ex) = c and B' x (ex) = d . Suppose first that fj, Ct( j(X) > /i c / jC (X). Then n c (X) > 
n c i{X). We define a coupling ^x of ttb x and ttb 1 as follows. Let (C,C) be a pair 
of colourings drawn from such that C is drawn from ttb x and C from 7r^/ . We 
have Pr[C(v x ) = c] = 0, Pr[C>x) = c'] = 0, Pr[C> x ) = c] > Pr[C(v x ) = d] and 
Pr[C"(wx) = i] < Pr[C(v x ) = i] for % e Q \ {c,d}. We pair up colourings in (C,C) 
such that C(v x ) = C'{v x ) when C'(v x ) = % for i e Q \ {c,d}. Then C{v x ) ^ C'{v x ) 
only when C'(v x ) = c. Thus, Pr[C(w x ) ^ C"(u^)] = At c ,c'(^) and < Vc,A x )- 

Suppose second that /i c ' C (X) > /i cc '(X). Similarly to above, v{X) < /i c / c (X). Thus, 
u(X)<ii(X). 

□ 

Suppose X is an edge-boundary pair and c = Bx{ex) and d = B' x (ex)- In order 
to obtain sufficiently good upper bounds on v{X) we use the previous lemma together 
with Lemma 12 below, which we first describe in words. Suppose we want to upper- 
bound fj, Ci( j(X). The idea is to pick a subregion R' C R x that contains vertex vx- Then 
we compute the maximum value of fi CyC > for that subregion, where we maximise over 
colourings of the boundary of R' that are identical to Bx on the overlapping boundary 
edges SRx H £R '. This maximum value is an upper bound on /i C;C i(X). Note that 
Goldberg, Martin and Paterson [12, Lemma 13] gave a similar lemma in terms of fi(X). 
However, in this paper it is crucial to be precise about the order of the colours c and d 
in Hc, c >{X). 

Lemma 12. Suppose that X is an edge-boundary pair and let c = B x {ex), c ' = B' x (e x )- 
Let R' be any subset of Rx which includes vx- Let S be the set of edge-boundary pairs 
X' such that Rx< = R' , the distinguished edge ex 1 = ex, and for the boundary colourings 
Bx> and B' x , we have Bx>(e) = Bx{e) and B' x ,(e) = B' x (e) on e e SRx H £R' . Then 
^c,d{X) < maxx'es Hc,c>(X'). 

Proof. Let X be an edge-boundary pair and let c = Bx{ex) and d = B' x (ex)- For a 
subregion R' C R x that contains Vx, let H = Rx\R'- For i e Q\{c} and 9 e £Ih, let 
riifi denote the number of colourings in Q Rx (B x ) which colour vx with colour i and H 
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(a) 




(b) 




(c) 




Figure 3: (a) An extended region TZ. A non-shaded vertex is labelled "in" , and a shaded 
vertex is labelled "out", (b) A region R with a marked edge e G E(R). (c) We see that 
TZ matches R with respect to edge e in R. 



with colouring 9. For 9 G Jiff, let n Ci denote the number of colourings in Q Rx (B' x ) which 
colour vx with colour c and H with colouring 9. Let N CtC > t g = J2 ie Q\( cc /) n i,e- Then 



n c (X) 



n c (X) + N CtC ,(X) Eeen H (nc,e + N c , dfi ) 



< max 



een H n c9 + N c 



< max/i C;C ,(X'). 



x'es 



To see the last inequality, take any 9 G Qh and construct the edge-boundary pair X' in 
S with the following parameters: Rx> = R', Bx> = Bx on £Rx H fit!' and -B^-, = B' x on 
H £ R' . For each boundary edge e G £ i?' such that e ^ £-Rx, let B x >(e) = B' x ,(e) = 
9(v), where vertex v G H is the endpoint of e in <9-R'. Now, 



7l f: 



n c (X') 



□ 



5.1 Extended regions 

It will be convenient to introduce the notion of an extended region 71, which is a region 
with the following additional information: (i) Every vertex in TZ is labelled either "in" or 
"out" , and (ii) one of the boundary edges of TZ is referred to as the designated edge. 

An extended region TZ and a region R are matching with respect to an edge e G E(R) 
if there is a way of overlapping R with TZ such that the designated edge of TZ coincides 
with the edge e, and every vertex that is labelled "in" in TZ coincides with a vertex 
that is in R, and every vertex that is labelled "out" in TZ coincides with a vertex that 
is not in R. When illustrating extended regions in the figures, we let non-shaded faces 
represent vertices that are labelled "in", and we let shaded faces represent vertices that 
are labelled "out" . We mark the designated boundary edge with a short and thick line 
segment. Figure 3 illustrates how an extended region TZ matches a region R with respect 
to an edge e. Note that the overlapping takes place under any rotation or reflection of 
the regions. 

Suppose TZ is an extended region. An extended region TZ' is an extended subregion of 
TZ if TZ' is obtained from TZ by removing vertices, except for the vertex that is incident 
to the designated edge. The labelling of the vertices in TZ' is identical to the labelling of 
the same vertices in TZ. 
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(a) 



(b) 




(c) 



Figure 4: (a) The extended region 1Zm (1 2) • (b) The extended region 7Zm (3 4) • (c) Labelling 
of vertices and edges. 






Figure 5: (a) The extended region Tip. (a)-(c) Three of the extended regions 
71 Fx, ■ ■ ■ ,71f 4720 - All vertices are labelled "in". 



5.2 A collection T of edge-boundary pairs 

Let 7Zm (1 2) be the extended region in Figure 4(a) and let 7Zm {3 4) be the extended region 
in Figure 4(b). Let M( 12 ) be the set of edge-boundary pairs X such that Rx and TZm^ 2) 
are matching with respect to ex- Let M 4 be the set of edge-boundary pairs X such that 
R x and 7^m (34) are matching with respect to ex- Let X be an edge-boundary pair and 
suppose c = B x (ex) and d = B' x (e x )- Let v be the vertex that is a neighbour to both vx 
and wx, let e be the edge between wx and v, and let d be the edge between v and vx (see 
Figure 4(c)). The three sets M x C M^ 2 ), M 2 C M (li2 ) and M 3 C M 4 of edge-boundary 
pairs are defined as follows. 

• X G Mi if f G -Rx and either // C;C / > /v, c and Bx(d) = c, or /v, c > // C)C / and 
B x (e) = d. 

• M 2 = M ( i )2 )\Mi. 

• X G M 3 if t> ^ i?x and either S x (e') = c or S x (e') = c'. 

Let Tip be the extended region in Figure 5(a). For / G {1, . . . ,4720} we define the 
extended region lZ Ff such that it is an extended subregion of Tip. Note that all vertices in 
1Zp ; are labelled "in". The remark on page 16 explains why we define exactly these 4720 
extended regions. Due to the large number of extended regions we only illustrate three 
of them here (Figure 5). For / G {1, . . . , 4720}, let Ff be the set of edge-boundary pairs 
X such that Rx and lZp f are matching with respect to edge ex- For m G {1, . . . ,4}, 
let Ff m = Fj fl M m . Let T be the collection of all sets F^ m . One of the extended 
regions TZp 1 , . . . ,TZp 4720 is defined to contain only the single vertex that is incident to 
the designated edge. Hence any edge-boundary pair X is guaranteed to belong to at 
least one of the sets in T . Note that many of the sets F^ m are empty. For instance, if 
TZp i is the extended region in Figure 5(b) for some % G {1, . . . ,4720} then obviously no 
edge-boundary pair X can belong to both Fj and M 4 . Hence F i>4: = 0. 
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5.3 The constants fiF fm 

For / G {1, . . . , 4720} we define F' f to be the set of edge-boundary pairs X such that the 
vertices of Rx are exactly those oITZf^ ex is the designated edge of 1Zf } , B x {ex) = 1, 
B' x (ex) = 2, and the number of colours q — 5. For an edge-boundary pair X G F'j, let v 
be the vertex that is a neighbour to both vx and let e be the edge between wx and 
f , and let e' be the edge between v and fx (see Figure 4(c)). Suppose first that 1Zm (1 2) is 
an extended subregion of Tl Ff . Then we define Fj x C F£ to be the set of edge-boundary 
pairs X E Ff such that B x {e) = 1, we define F£ 2 C Fj- to be the set of edge-boundary 
pairs X <E F'j- such that Bx(e) = 2, and we define Fj 3 = F'j A = 0. Suppose second that 
1Zm, 12 ) 1S n °t an extended subregion of TZp f - Then we define F'^ — F'^ 2 — 0, we define 
^ -F/ to be the set of edge-boundary pairs X G Fj- such that either Bx(e') = 1 or 
Bx(e') = 2, and we define F' f A = F' f . Now, for / e {1, . . . , 4720} and m G {1, . . . , 4}, we 
define 

A*F /im = max ^1,2 (X) 

A t-T f 

if F' Lm ± 0, and ^ /im = if F' f m = 0. 

Lemma 13. Suppose q = 5, f G {1, . . . , 4720}, m G {1,...,4} and F f , m ^ 0. Then 
v(X) < fiF fm for every edge-boundary pair X G F^ m . 

Proof. Suppose / G {1, . . . ,4720} and m G {1, ... ,4} such that F fiTn ^ 0. Let X be 
an edge-boundary pair in F/, m . Let f be the vertex that is a neighbour to both vx and 
wx, let e be the edge between wx and v, and let e' be the edge between v and vx (see 
Figure 4(c)). From Lemma 11 we have that v(X) < fi(X). In order to upper-bound 
li{X) we may assume without loss of generality that Bx{ex) = 1 and B' x (ex) = 2. 

Suppose first that m — 1. Without loss of generality we may assume that B x {e) = 1 
and hence / Ui j2 (X) > /i 2 ,i(X). Then /i(X) = fj, lt2 (X). Let R' be the subset of Rx such 
that the vertices of R' are exactly those of lZp f . Let S be the set of edge-boundary pairs 
X' such that Rx> = R', the distinguished edge ex> = ex, and for the boundary colourings 
B x > and B' x , we have B x >{e") = B x (e") and B' x ,(e") = B' x (e") on e" G £R x n£R'. Note 
that S C We have 

/ii, 2 (X) < max/ii i2 (X') < max ^i, 2 (X') = fx Ff m , 

where the first inequality is from Lemma 12. 

Suppose second that m — 2. Without loss of generality we may assume that B x (e) = 2 
and hence fj,i j2 (X) > fi 2 ,i(X). Proceeding as above we see that /jl(X) < [±F fm - Now 
suppose m — 3. Without loss of generality we may assume that B x (e') = 1 or B x (e') = 2 
and ni )2 (X) > fi 2> i(X). Proceeding as above we see that fi(X) < f^F fm - Lastly, for m = 4 
we make no assumption on the colour of edge e' and again we see that fi(X) < Hf } m - 

□ 

5.4 A collection A of edge-boundary pairs 

Let IZa be the extended region in Figure 6(a). For a G {1, . . . , 342} we define the extended 
region 1Za u to be a subregion of TZa- The extended regions are defined such that for any 
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Figure 6: (a) The extended region TZa- (a)-(c) Three of the extended regions 
TZai 1 ■ ■ ■ ) TZa 3 42 • 

edge-boundary pair X, the region Rx matches exactly one of 7Za ± , ■ ■ ■ , TZa 342 with respect 
to edge ex- The remark on page 16 explains why we define exactly these 342 extended 
regions. In Figure 6 we illustrate three of the 342 extended regions. For a G {1, . . . , 342}, 
let A a be the set of edge-boundary pairs X such that Rx matches TZA a with respect to 
edge ex- Furthermore, for m G {1, . . . , 4} we define A a ^ m = A a n M m , and define A to be 
the collection of all sets A a ^ m . Note that many of the sets A a ^ m are empty. 

5.5 Exponential decay 

AsetiSC.AxjFx^lx^lx^lis called an (A, J 7 ) -set if the following is true about S: For 
every set A a , m G A, every edge-boundary pair X G A a , rn , and every two distinct colours 
c,d e Q such that p™ n (c, c') > 0, there is a 5-tuple (A atTn , F fim , A aumi , A a2>m2 , A a3jm3 ) in 
S, such that X G i 7 /,™., and for « G {1, 2, 3} the edge-boundary pair Xj(c, c') constructed 
recursively in the tree T x belongs to A a . >m .. For values of % such that X^c, d) = 0, 

A ail rrii 0- 

Suppose e G (0,1) is a constant. An (A,J-)-set S is good with respect to e if the 
following is true: For i G {1, . . . , 342} and j G {1, . . . , 4} there is a constant a:^ such 
that otAij > 1/(1 — e) if 7^ and a:^.., > if = 0, and for every 5-tuple 
(A at m, Ff ;7n , A ai>mi , A a2jrrl2 , A a3tTn3 ) in S, 

VF f , m (aA ai , mi + Oi Aa2>m2 + a A a . i>mz ) < «A , m (l - e). (3) 

Lemma 14. Suppose q = 5, e G (0, 1) a constant, and S is an (A,J r )-set that is good 
with respect to e. Then there is a constant a > such that T d (X) < a(l — e) d for all 
edge-boundary pairs X . 

Proof. Since S is good with respect to e, there are constants ctA a , m , a G {1, . . . , 342} and 
m G {1, . . . , 4}, such that Equation (3) is satisfied for every 5-tuple in S. For A a>m G ^4, 
let T d (A a>m ) denote the maximum of T d (X) over all X G A a>m . Remember r d (0) = for 
d > 1. In order to show that there is a constant a such that T d (X) < a(l — e) d for every 
edge-boundary pair X, we will show that T d (A a ^ m ) < a J 4 a m (l — e) d for every non-empty 
set A a ^ m G A. Then we let a be the maximum of ctA a . m over all a G {1, . . . , 342} and 
m G {1, . . . , 4}. Note that any edge-boundary pair X belongs to at least one of the sets 
in A. 

Consider any non-empty set A a>m G A and any edge-boundary pair X G A a ^ m . We 
are going to show that T d (X) < a J 4 a , m (l — d) by induction on d. We start with the base 
case d—1. Since a Aa ^ m > 1/(1 — d), we have 

Y\X) = u(X) < n(X) < 1 < a Aa , m {l - e), 



14 



where the first inequality is from Lemma 11. Now consider the inductive step. We repeat 
Equation (2): 

T d (X)= J2 PifMjt^M). ( 4 ) 
c,c'eQ, i=i 

where Xj(c, c') is the edge-boundary pair constructed recursively in the tree Tx- Here 
Q = {1, . . . , 5}. For every two distinct colours c,c' G Q such that p™ in (c, c') > 0, we know 
that there is a 5-tuple (A ajTn , F f , m , A aumi , A a2jrjl2 , A a ^ m3 ) in S such that X G F f , m and 
Xi(c, d) G A aumi , where i G {1,2,3}. If the i-th neighbour of Vx is not in Rx then we 
have A aumi = 0. By the induction hypothesis we have 

X>^(X,(c,c')) < j^T d -\A a ^ nH ) < j^a Aw {l-e) d -\ (5) 

i=l i=l i=l 

Using Equation (4) with Equation (5) gives 

3 3 

c,c'eQ, i=l i=l 

3 

< J>a w (1 - e)^ 1 < aAaim (l - e) d , 
i=i 

where f(-X") < //f. is from Lemma 13, and the last inequality follows from Equation (3). 

□ 

The next lemma is proved by computer assistance and we will explain the details in 
Section 8. 

Lemma 15. Suppose q = and e = 1/1000. Then T d {X) < 5(1 - e) d for every edge- 
boundary pair X . 

Proof. In order to prove this lemma we use computer assistance. The computerised steps 
are to first calculate all the constants //F /m , then generate an (A, J 7 ) -set S that is good 
with respect to e — 1/1000. This last step is broken into the following steps. First we 
generate an (A, J 7 ) -set S. Then, for every 5-tuple (^4a,mj Ff,mi A aitmi , A a2 ^ m2 , A a3jJ7l3 ) in 
S, we add the inequality 

Up, (a A +&A + &A ) < CtA (1 — S) 

to a linear program. The unknowns in this linear program are the variables «A am - A 
solution to the linear program is found with a.A a . m G [2, 5] for A a , m ^ and aA a>m G [0, 5] 
for A a , m = 0. Hence S is good with respect to e. By Lemma 14 it follows that T d (X) < 
a(l — e) d for every edge-boundary pair X, where a > is a constant. From the proof of 
Lemma 14 we see that we can choose a to be the maximum of all a^^, which is 5. □ 
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Remark. One probably asks why the sets in A and T are the sets we use to prove 
mixing. The sets in A and or the extended regions TZ^ and 7Zp l to be more precise, 
have arisen from a lengthy process of trial and error and experiments. One part of the 
proof of Lemma 15 above is to find a solution to a linear program. If the values (J>F fm 
are too large then there will be no solution to this linear program. In order to obtain 
smaller values ^F fm we must increase the size of the regions 1Z Fi . Small extended regions 
IZa, contain only little information about which vertices are in and not in the region Rx 
for an edge-boundary pair X G A,. In particular, with small regions IZ^. we quickly lose 
information about which vertices are in and not in the regions Rxi(c,c') f° r the recursively 
constructed edge-boundary pairs X i (c,c / ). Thus, too small extended regions 1Z a % will 
result in a linear program that is too small and has no solution. We started with a 
few small extended regions TZai and TZp t and slowly increased the sizes of them until we 
obtained a linear program that could be successfully solved. We let the regions grow in 
a way that seemed reasonable based on experiments and intuition. 

6 Strong spatial mixing 

Lemma 16. Suppose e = 1/1000 and q = 5. Suppose X is a vertex-boundary pair and 
R' Q Rx- Then there is a coupling ^ x of t^b x an d n B' x suc h that 

veR' ^ ' 

Proof. First suppose that wx has a neighbour y ^ Rx- Let k — \E\ < 3, where E = 
{ei, . . . , Ck\ Q £Rx is the set of boundary edges incident to wx- Label the edges in E 
clockwise around wx so that edge {wx, y) appears between edge and e± when traversing 
edges around wx in clockwise direction. This guarantees that and ej are adjacent only 
if % and j differ by 1. 

For i — 1, . . . , k, let be the edge-boundary pair consisting of region R Xi = Rx, the 
distinguished edge ex t = ei, and boundary colourings Bx t and B' x ■ For every boundary 
edge e = (w,v) € SRx t \E, where w G dRx, we have BxX e ) — B' Xi (e) = Bx(w). The 
colours of the edges in E are assigned as follows. 

. B Xi (ej) = B'xXej) = B' x {w x ) for j = 1, . . . , % - 1, 

• B x .(ej) = Bx(w x ) and B' x .(ej) = B' x {w x ) for j = i, and 

• B Xi (ej) = B'x.iej) = Bx{w x ) for j = i + 1, . . . , k. 

By Lemma 10 there is a coupling of 71b x an d tt_b' such that 

veR' d>d(e x .,R') 

Let ^ x be the coupling of txb x and hb' x defined by composing the couplings ^i, . . . , ^k- 
More precisely, in order to choose a pair (oojOfc) of colourings from ^x, first draw the 
pair (oo,o"i) from ^i. Say ctq = x and o\ = x\. Then choose the pair (ai,^) from the 



16 



conditional distribution ty 2 , conditioned on o\ = %\. Say a 2 = x 2 . Then choose the pair 
((72,(73) from the conditional distribution ^3, conditioned on a 2 = x 2 , and so on. Hence, 
<t is drawn from n Bx = n Bx and is drawn from it b i = n B i . By the construction of 
the coupling it follows that if the colour of a vertex v G Rx differs in a pair (cr , a^) 
drawn from ^x then it must differ in at least one of the pairs ((7j_i, o"j) drawn from 
where % = 1, . . . , k. Using Equation (6) and Lemma 15 we have 

k k 

E nu XtV ] < E E E ^ = EE E [wi 

vGi?' veR' i=l i=l veR' 

<t £ r"«)<£ £ 5(1-^ 

i=l d>d{e X -,R') »=1 d>d(w x ,R>) 




Now suppose all neighbours of w# are in i?^. Breaking the discrepancy at vertex 
into edge-boundary pairs Xj as above is not possible because the induced edge-boundary 
pairs are not valid with respect to the colouring of adjacent boundary edges. 

Let u G Rx be a neighbour of wx- Suppose u R'. Let Rx, u = Rx\{u} be the 
region Rx after removing vertex u. For c G Q, let Bx, c be the colouring of the vertex- 
boundary dRx, u such that for all v G dRx H dRx, u , Bx,c{v) = Bx(v), and Bx, c (u) = c 
(if ii G dRx, u )- Similarly, for c' G Q, let £>^, c , be the colouring of the vertex-boundary 
dRx, u such that for all v G dRx H 8Rx, u , B'x c'( v ) = B'x( v )i an d ^ C '( M ) = c '- Note that 
the colourings and £>^ c / can differ on up to two vertices, namely on vertex wx and 
u. We break the difference in the (up to) two vertices wx and u on the boundary dRx, u 
into differences in the edges that bound them. 

Let k = \E\ < 6, where E = {e±, . . . ,ek} C £Rx, u is the set of boundary edges 
incident to Wx or u. Label the edges in E clockwise around wx and u so that e& and e\ 
are not adjacent. Such a labelling is always possible since wx and u are neighbours. This 
guarantees that and tj are only adjacent if i and j differ by 1. 

Let c E Q and c' G Q be two (not necessarily different) colours. Similarly to above, 
for i — 1, . . . , k, let Xi be the edge-boundary pair consisting of region Rx t = Rx,u, the 
distinguished edge ex { = and boundary colourings Bx t and B' x .. The colourings Bx { 
and B' x . are defined similarly to above, as a sequence of colourings differing only on the 
distinguished edge e^. That is, for a boundary edge e = (w, v) G £Rx, u , where w G dRx, u , 
we have Bxi(e) = £>a\c(w) and B' x (e) = <B^ c ,(u>). Let Vf^ be a coupling of 7Tb x and 7rg/ 
such that Equation (6) is satisfied, which possible due to Lemma 10. We now construct 
a coupling ^ x of ~kb x and 7Tg^, in the following way. 

Let 'fy'x be any coupling of 7Tb x and t^b' x - Let (C, C") be the random variable corre- 
sponding to the pair of colourings drawn from \I> x (yet to be constructed). We will choose 
the colour of u in C and C according to ^' x . Let c and d be the colour of u drawn from 
W x . Let ^x,c,c' be a coupling of n Bx c and 7r^ r To complete the construction of ^> x 
we colour the remaining vertices in Rx by choosing two colourings from ^x,c,&- The 
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coupling ^x, c ,c' is constructed by composing the k couplings as above. We have 

k k 

E e[i**,«] < E E = E E E [ 1 **.«] 

uGi?' vGR' i=l i=l uefl' 

<£ £ r"(A-,)<£ J] 5(1 -e)' 

i=l d>d{e Xv R') »=1 d^dCwAr.-R')- 1 

= E -(! - e)^'^" 1 < -7^(1 - e) 



30 • ^d(w x ,R') 



i=l 



where the —1 in u d(wx, R') — 1" comes from the fact that the distance from the discrep- 
ancy edge to R' may be one less than d(wx,R')- Since we sum over all distances 
greater than or equal to d(wx, R') — 1, and (1 — e)° = 1, we note that the bound also 
holds when u e R' . □ 

We now prove Theorem 4 of strong spatial mixing for q = 5 colours. 

Theorem (4, repeated). The system specified by proper 5-colourings of the kagome lattice 
has strong spatial mixing. 

Proof. Consider the vertex-boundary pair X such that, from Definition 1 of strong spatial 
mixing, we have Rx = R, Bx = B, B' x = B' and wx = w. Let R' be any subregion of R. 
The total variation distance between hb,r' and Kb',r' is upper-bounded by the probability 
that R' differ under any coupling ^ of hb,r' and ttb',r'- This probability is upper-bounded 
by YlveR' ^[-'-*' v \- Using the coupling ^ x in Lemma 16, we have 

dTv(ir B ,R>^&,R') = d TY (7TB x ,R',7TB' x , R ') < E^ 1 ^^] < a\R'\(l - e)^^ , 

veR' 

where e = 1/1000 and a = 30/(e(l - e)). □ 



7 Rapid mixing 

The implication from strong spatial mixing to rapidly mixing Glauber dynamics is only 
known to hold for graphs of sub- exponential growth [25], meaning that the number of 
vertices at distance d from any vertex v is sub-exponential in d. This is an important 
property we make use of in the proof of rapid mixing in this section. For further discussion 
on this topic in general, see [12], in particular [12, Section 7.5]. 

Lemma 17. Let v e Vg be any vertex in the kagome lattice and let rid(v) denote the 
number of vertices at distance d from v. Then n d {y) G 0(d). 

Proof. Recall the definition of the kagome lattice in Section 1.1, in particular Figure 1. 
First assume that v G Kdd- In order to derive lower and upper bounds on n^(f), we 
assume without loss of generality that v = (1,1) is the vertex at re-coordinate 1 and 
y-coordinate 1. Fix any positive integer d. 
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We first derive a lower bound on n d ((l, 1))- For each odd value of y G {1, . . . , d}, 
let (x,y) be the vertex at distance d from (1, 1) that is reached with the following path: 
(1, 1), (2, 2), (3, 3), . . . , (y, y),(y + 2,y),(y + 4,y)..., (x, y). Note that vertex (y, y) G V odd , 
and from (y, y) we go as far as possible to the right. Also note that there is no path from 
(1,1) to (x,y) that is shorter than length d. Thus, there are at least [d/2\ vertices at 
distance d from (1, 1), and we have n d ((l, 1)) > [d/2\. 

When deriving an upper bound on n d ((l, 1)) we will use two claims: 

Claim 1. For any two vertices (x,y\ ow ) and (x, ywgh), where 1 < yi ow < ytigh, the 
distance between (1,1) and (x, y\ ow ) is strictly smaller than the distance between (1,1) 
and (x,yhi g h)- We prove the claim by considering two cases: 

Case (i). Assume that x is odd, and hence both (x, y\ ow ) and (x, yhigh) are in V^d- 
Consider a shortest path from (1, 1) to (x, yugh)- The path must use a vertex (x pass , yi ow ) G 
Kdd at y-coordinate y low . From (x pass , y low ) we can reach (x,y iow ) in exactly \x - x pass |/2 
steps. The number of steps required to reach (x,yugh) from (x pass , y\ ow ) is strictly greater 
than \x — :r pass |/2 since some steps must be used to increase the y-coordinate so it will 
eventually reach yugt, and for each such up-move the ^-coordinate is increased/decreased 
only by 1. Thus, if x is odd then the distance between (1,1) and (x, y\ ow ) is strictly smaller 
than the distance between {x,yugh)- 

Case (ii). Assume that x is even, and hence both (x, y\ ow ) and (x, i/high) are in V cven . 
We will use the same argument as for odd values of x, only with the difference that we 
consider a vertex (x pSuSS ,yi ow — 1) G V a dd on a shortest path from (1, 1) to (x, t/high)- From 
{xpass, l/iow — 1) we can reach (x, y\ ow ) in at most [\x — rc pass |/2j + 1 steps, where the +1 
comes from the fact that we need to go up one ^-coordinate. The number of steps required 
to reach (x, ytigh) from (x pass ,yi ow — 1) is strictly greater than [\x — x pass \/2\ + 1 since 
some steps must be used to increase the ^-coordinate so it will eventually reach yugh, and 
for each such up-move the ^-coordinate is increased/decreased only by 1. Thus, also for 
even values of x we have that the distance between (1, 1) and (x, y\ ow ) is strictly smaller 
than the distance between (1, 1) and (x, yugh)- 

Claim 2. For any two vertices (x, y\ ow ) and (x, ytigh), where y\ ow < i/high < 1, the 
distance between (1, 1) and (x, yhigt) is strictly smaller than the distance between (1, 1) 
and (x, |/iow)- We prove the claim by using exactly the same reasoning as for Claim 1. 

Using Claim 1 and 2 we conclude that there are at most two vertices (x, y) and 
(x,y'), with the same x-coordinate, at distance d from (1, 1). The leftmost vertex that 
is at distance d from from (1,1) is (1 — 2d, 1). It is reached by making d consecutive 
left-moves. Similarly, the rightmost vertex at distance d from (1, 1) is (1 + 2d, 1). Thus, 
the ^-coordinate of any vertex at distance d from (1,1) is in the set {1 — 2d, . . . , 1 + 2c?}, 
and hence there are at most 2 x (4d +1) = 8d + 2 vertices at distance d from (1, 1). 
That is, n d ((l, 1)) < 8d + 2. We have now showed that for any vertex v G Kdd, [d/2\ < 
n d (v) <8d + 2. 

It remains to derive upper and lower bounds on nd(v ) for v G V even . Without loss of 
generality we assume that v — (0, 0) is the vertex at re-coordinate and y-coordinate 0. 
Fix any positive integer d. 

We derive a lower bound on n^((0, 0)) in the same way as when v = (1, 1). For each 
odd value of y G {1, . . . , d}, let (x, y) be the vertex at distance d from (0, 0) that is reached 
with the following path: (0, 0), (1, 1), (2, 2), ... , (y, y), (y+2, y), (y + A, y)...,(x,y). Thus, 
there are at least [d/2\ vertices at distance d from (0, 0), and we have nd((0, 0)) > [d/2\. 
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We now derive an upper bound on n d ((0, 0)). Vertex (0, 0) has exactly four neighbours: 
(1, 1), (1, —1), (—1, —1) and (—1, 1), which are all in V^d- The shortest path from (0,0) 
to any vertex at distance d from (0, 0) must use one of these four vertices. Thus, an upper 
bound on the number of vertices at distance d from (0,0) is n<i((0, 0)) < n d -i((l, 1)) + 
rid_i((l, —1)) + 7i(j-i((— 1, —1)) + 7i(j_i((— 1, 1)). From the upper bound above we have 
that there are at most 8(d — 1) + 2 vertices at distance d — 1 from a vertex in Kdd- Hence 
there are at most than 4 x (8(d — 1) + 2) = 32d — 24 vertices at distance d from (0, 0), 
and we have n d ((0, 0)) < 32d - 24. 

Finally, for any vertex v G Vg and any positive integer d we have shown that [d/2\ < 
n d (v) < 32d - 24. □ 

For a vertex v G Vg and an integer d > 0, let Ballot") denote the set of vertices that 
are at most distance d from v. Thus we have Ballo(n) = {v }. 

Lemma 18. For any rea/ number a > inere is an integer d > sfica £/ia£ 

|9Ball d (n)| 



|Ball,(n)| -°' 



uniformly in n G Vg. 



Proof. Let n be a vertex in Vg and let a > be a real number. For an integer d > 0, let 
n<j(n) denote the number of vertices at distance d from n. By Lemma 17, n^n) G O(d). 
We have |<9BalLj(n)| = n d+ i(v) G 6(d) and |BalLj(v)| = Ylt=o n d{ v ) e @(d 2 ). Hence there 
is an integer d > such that |9Ball d (n)|/|Ball d (n)| < a for d > d . □ 



7.1 The Markov chain Md 

In order to analyse the mixing time of the Glauber dynamics we first define a similar 
Markov chain that corresponds to heat-bath dynamics on small subregions instead of 
single vertices. For a region R, vertex v G Vg and integer d > 0, let R% — R n BalLj(n). 
Let R d = {v G Vg | R% 7^ 0}. For a region i?, go-colouring i3 of OR and integer d > 0, 
we define the heat-bath Markov chain M. d as follows. The state space is VLr{&) and a 
transition from a state a is made in the following way: First choose a vertex v uniformly 
at random from R d . Let B d be the colouring of dR d induced by a and £>. To make 
the transition from a, recolour the vertices in Rf, by sampling a colouring from 7Tgd, the 
uniform distribution on proper colourings of the region R d that agree with dR d . As for 
the Glauber dynamics, the stationary distribution of M. d is vr^. Since Ballo(n) = {n}, 
Glauber dynamics is M.q. In order to prove rapid mixing of the Glauber dynamics, we 
will use the mixing time of Ai d for some constant d and use a Markov chain comparison 
method to infer rapid mixing of Aio- 

To establish the mixing time of A4 d we use path coupling, due to Bubley and Dyer [3]. 
Let o~\ and o~i be two states of A4 d , where d is to be specified. Using the path-coupling 
method, we only need to consider two colourings U\ and a 2 that differ on exactly one 
vertex, which we refer to as w. That is, the Hamming distance between a x and a 2 is 1. 
Let Ai d make a transition from o~\ to o~[, and from o~ 2 to a' 2 . We want to correlate 
(or couple) these two transitions such that the expected Hamming distance between a[ 
and a' 2 is less than 1. If we can do this then we use the path-coupling theorem (see 
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for instance [3, 7]]) to infer the mixing time of Aid- It is possible to construct such a 
coupling of the transitions provided d is sufficiently large. The idea is that we update 
the same vertices R% in both the transition from o~\ to a[ and o~ 2 to a' 2 . If the vertices 
we update do not include w, and w is not in dR%, then we choose the same colouring of 
R% in both transitions, and hence the Hamming distance between a[ and a' 2 remains 1. 
If the vertices R% we update contain w then again we choose the same colouring of R% 
in both transitions, and the Hamming distance drops to 0. The only situation when the 
Hamming distance can increase is when w is on the boundary dR„ of the vertices R% we 
update. In this case we use the coupling in Lemma 16 to colour the vertices in R%. This 
guarantees that the expected Hamming distance between a[ and a' 2 will only increase 
by at most a constant K = 30/(e(l — e)). Due to Lemma 18 we can choose a radius d 
such that the ratio of the probability of having w G dR% and the probability of having 
w G is arbitrarily small. Thus, we choose d such that the probability of decreasing 
the Hamming distance by 1 is so much bigger than the probability of increasing it by K 
that the expected Hamming distance between a[ and a 2 is less than 1. The exact details 
of how to achieve this is explained in Sections 7.1 and 7.2 in [12]. In Section 7.2 in [12] 
a proof of the following lemma is found. Note that the notation in [12] differ slightly 
and of course we make use of Lemmas 16 and 18 as explained above rather than using 
equivalent lemmas in [12]. 

Lemma 19. Suppose q = 5. There is an integer d > such that the Markov chain Aid 
is rapidly mixing on any region R under any q -colouring B of dR. The mixing time 
T M d {$) £ 0{n\ogj), where n is the number of vertices in R. 

7.2 Rapidly mixing Glauber dynamics 

We will compare the mixing time of the Markov chain Aid and the Glauber dynamics 
Aio by using a method of Diaconis and Saloff-Coste [4]. Their method has been used 
before by Goldberg, Martin and Paterson in [12] to compare the mixing time of Aid and 
Aio under the assumption that q > A + 2, where A is the maximum degree of the lattice. 
Here we consider q = 5 on the kagome lattice (A = 4) and therefore we cannot make 
direct use of the comparison in [12]. Next we review the comparison described in [12] and 
provide a proof of rapidly mixing Glauber dynamics with q = 5 colours. For a survey on 
Markov chain comparison in general, see [6]. 

Let Pd and Pq denote the transition matrix for the chain Aid and Aio, respectively. 
For % G {0, d}, let be the set of pairs of distinct colourings (a±, o~ 2 ) with Pi(cri, o~ 2 ) > 0. 
The set Ei can be thought of as containing the edges of the transition graph of Aii, and 
hence we sometimes refer to a pair in Ei as an edge. For every edge (o"i,o"2) G Ed, let 
V ai ,a 2 be the set of paths from U\ to a 2 using transitions of Ai . More formally, let V ai>er2 
be the set of paths 7 = (<7i = 9 , 9\, . . . , 9 k = a 2 ) such that 

(1) each (9i,0i+i) is in E , and 

(2) each edge in E appears at most once on 7. 

We write I7I to denote the length of path 7. So, for example, if 7 = (9 , ...,9k) we have 
I7I = k. Let V = ^{a 1 ,<T 2 )&E d V<T 1 ,<T 2 t> e the set of all paths for all edges in E d . 
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A flow is a function from V to the interval [0, 1] such that for every (a±, <r 2 ) G E d , 



J~] 0(7) = P d (fJi, (7 2 )7r B (o-i). 
For every (#1, # 2 ) G i?o, the congestion of edge (#1, #2) in the flow is the quantity 

v ; v ' ; 7eP:(ei,e 2 )€7 
The congestion of the flow is the quantity 

^(0) = ™ax A? l5 e 2 (0). 

(6»i,6» 2 )G-Bo 

Theorem 20 below describes how the mixing times of Aid and .Mo are related. A proof 
of this theorem can be found in [6, Observation 13]. As pointed out in [12], this theorem is 
similar to Proposition 4 of Randall and Tetali [19] except that [19, Proposition 4] requires 
the eigenvalues of transition matrices to be non-negative. Both results are based closely 
on the ideas of Aldous [2], Diaconis and Stroock [5], and Sinclair [21]. Let TM d {$) be the 
mixing time of and let t Mq (5) be the mixing time of the Glauber dynamics .Mo- 
Theorem 20. Suppose that is a flow. Let p = min ee Q R ^ P (9, 9) and assume that 
p > 0. Then for any < 5' < \ 



tm (S) < ln -r— max 

■ 7T m i n 



where 7r min = min (T€ n J? (B) ir B {o). 



Lemma 21. Suppose that there is a flow such that the congestion A(<p) e 0(1). Then 
the mixing time of the Glauber dynamics _M on a region R is r Mo (5) G 0{n{n + log |)) ; 
where n is the number of vertices in R. 

Proof. From Definition 2 of Glauber dynamics, p = ming e ^ R ^ P (6,6) > 1/q. Suppose 
5' = 1/n. Then by Lemma 19 we have T>| d (<5') G O(nlogn). With A((j>) G 0(1), 
Theorem 20 gives 

T Mo (5) < In — !— • 0(1) • 0(n) = 0(n(n + log -) 

■ TTmin 

since 7r min > l/q n and hence ln(l/7r min ) G 0(n). □ 

In order to establish the mixing time of the Glauber dynamics Aio by applying 
Lemma 21 we have to construct a flow such that the congestion A{<p) G 0(1). Given a 
(/-colouring a of a region R and a go-colouring B of dR, a single-vertex update of a vertex 
f G i? is a recolouring of v to a colour c E Q such that no neighbour of v has colour c 
in either a or £>. Suppose R is a region and <Ji and o" 2 are two proper 5-colourings of i? 
that differ on m vertices. The next two lemmas tell us how a series of 0(m) single- vertex 
updates applied to <Ji can transform G\ to cr 2 - This sequence of single- vertex updates will 
be used when constructing the flow 0. 
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Figure 7: In every proper 5-colouring of the region in (a) there is a vertex that has two 
neighbours with the same colour. 

Lemma 22. Consider the region in Figure 7(a). In every proper 5-colouring of this 
region there is a vertex that has two neighbours with the same colour. 

Proof. Suppose a is a proper 5-colouring of the region in Figure 7(a) such that no two 
neighbours of a vertex in the region have the same colour. We will show that this leads to 
contradiction. Without loss of generality we may assume that five of the vertices have the 
colours specified in Figure 7(b). A vertex is labelled with its colour. It follows that the 
two vertices adjacent to the vertex coloured 5 must have colour 3 and 4, otherwise there 
would be a vertex that has two neighbours with the same colour. Similarly, the vertices 
adjacent to the vertex coloured 3 must have colour 1 and 5, and therefore the two bottom 
left vertices must have colour 2 and 4 in a. Figure 7(c) illustrates this fact, where a square 
contains the two colours of the two vertices it is overlapping. From the two left squares 
we see that the colour 4 must be on the vertices that are as far apart as possible. Thus, 
a must agree with the colouring in Figure 7(d). Figure 7(e) illustrates how other vertices 
of the region must be coloured in a, and Figure 7(f) shows the necessary colouring of the 
four rightmost vertices at the top. To finish the proof we note that it is impossible to 
assign colours to the two leftmost vertices at the top without introducing a vertex such 
that two of its neighbours receive the same colour. □ 

Lemma 23. Let R be a region of the kagome lattice and let B be the 0-colouring of the 
boundary OR. Suppose that q = 5 and let o\ and a 2 be any two proper q-colourings of 
R that differ on m vertices. We can go from a x to a 2 by applying a series of 0{m) 
single-vertex updates. 

Proof. Let v G R be a vertex on which o\ and a 2 differ. We will show how to recolour v 
to the colour it has in <t 2 by doing at most a constant number of single-vertex updates. 
A vertex in R that has the same colour in both o\ and a 2 will not change colour after 
v has been updated. First we analyse situations where no boundary vertices in dR are 
involved. We note at the end of the proof that if boundary vertices are present, then it 
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Figure 8: The colours 1 and 2 are going to swap place. Lemma 23 guarantees that this 
can be done with a constant number of single-vertex updates. 

only makes it easier to recolour v. That is, assume for now that all vertices we consider 
belong to the region R. The proof goes through a series of cases. 

If possible, simply recolour v to the colour it has in cr 2 . If this is not possible then 
there must be one or two neighbours of v that have colour o^f) i n It cannot be more 
than two such neighbours since U\ is a proper colouring. 

Without loss of generality, assume that o\{v) = 1 and a 2 (v) =2. If v has two 
neighbours with colour 2 in <j\ then we will first recolour one of these two neighbours to 
some other colour than 2. Let w be the neighbour of v with colour 2 that we are going 
to recolour. Note that U2(w) ^ 2 since 02 is a proper colouring. If possible, recolour w 
to some other colour than 2. If this is not possible then w is "locked" and must have 
three neighbours coloured 3, 4 and 5, respectively. In this case, first recolour v (which is 
possible since v has two neighbours with colour 2) and then recolour w to colour 1. Now 
only one neighbour of v has colour 2. We deal with this case next. 

Without loss of generality, assume that o"i(t>) = 1 and o" 2 (f) = 2, and exactly one 
neighbour w of v has colour 2 in <j\. Note that oiiyS) 7^ 2 since 02 is a proper colouring. 
If possible, recolour w to something else than 2 and then recolour v to 2. If this is not 
possible then w is "locked" and must have four neighbours (including v) with colours 1, 3, 
4 and 5, respectively, in o\. Without loss of generality, consider the region in Figure 8(a), 
which is a subregion of R. Call this region R' . The vertices of R' are labelled with their 
colours in cr 1 . The vertex with colour 1 is v and the vertex with colour 2 is w. We 
assume without loss of generality that the two neighbours of v that are below v are the 
two neighbours with colour 3 and 4 in o\. 

Three of the vertices in R' are given the colours a, b and c, which are to be determined. 
Since w is "locked", the colours a, b and c is any permutation of the colours 3, 4 and 5. 
If a is 3 or 4 then we recolour v to 5 and then recolour w to 1, and then recolour v to 2. 
If this is not the case then a must be 5, and hence the colours b and c are 3 and 4 in any 
order. Figure 8(b) illustrates this. We now analyse this case. 

We will use Lemma 22 to show that we can recolour v to 2 without changing the 
colour of any other vertex except w (which will be recoloured to 1). Consider Figure 8(c) 
which illustrates the region R' extended with vertices in R. The vertices we extend R' 
with correspond to the region that we used in Lemma 22. From Lemma 22 we know that 
there must be at least one vertex u among the vertices we extend R' with such that u has 
two neighbours with the same colour. Let P be a shortest path from v to u such that the 
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path goes from v to the neighbour above that has colour 5 and then is entirely inside the 
region we added to R' . Figure 8(d) illustrates an example of such a path. The path is 
shaded in the figure. Suppose that the vertex u is chosen such that all vertices on the path 
P (except from u itself) are "locked" (have four neighbours of different colours). Note 
that if the vertex coloured 5 above v does not have four neighbours of different colours 
then we let u be this vertex and hence the path P consists only of the two vertices u and 
v. 

Suppose that the path P contains k vertices. Let c±, . . . , be the colours in o\ of the 
vertices from u to v along the path. That is, (J\{u) = c±, Ck-i = 5 and — Since u has 
two neighbours with the same colour, we recolour u from C\ to another colour c[. Now 
the vertex after u on P has two neighbours with the same colour (namely c'J, since all 
its neighbours had different colours before recolouring u. We recolour this vertex from c 2 
to c' 2 . We continue this recolouring procedure along the path P all the way to vertex v, 
which will be reco loured to 3. Note that the vertex above v which had previously colour 5 
now must have colour 3 or 4. We can now recolour w to 1 and then recolour v to 2. It 
remains to recolour the vertices on the path back to their original colours in o\. We do 
this by reversing the recolouring procedure, starting with the vertex above v, which is 
recoloured back to 5. When u is recoloured back to c\ we are done. 

We have now shown how a constant number of single-vertex updates are applied in 
order to recolour a vertex v to the colour it has in a 2 without changing the colour on 
vertices that have the same colour in a x and a 2 . 

We note that if any vertices involved in the recolouring procedure of v are boundary 
vertices then this will only make it easier. Note from the statement of the lemma that we 
assume that a boundary vertex has colour 0. As we have seen, the tricky situations arise 
when a vertex is "locked" with four neighbours of different colours (excluding colour 0). 
Such a vertex is tricky because we cannot just change its colour to another colour in 
Q = {1, . . . , 5}. A vertex that is adjacent to a boundary vertex can never be "locked" 
since there is always at least one colour in Q that it can be recoloured to. Thus, although 
the part of the proof above assumes that all vertices are in R, we note that the presence 
of boundary vertices only makes the recolouring procedure easier. Of course, depending 
on which vertex v we are going to recolour, and which neighbour w is "locked" , the path 
P might go in a direction that is different from the one in Figure 8(d). However, the 
same technique is applied in order to successfully recolour v. 

Finally, in order to transform U\ to a 2 , we recolour each vertex v at which U\ and a 2 
differ. For each such vertex it takes only a constant number of single-vertex updates to 
do so. Since u\ and a 2 differ only at m vertices, the total number or updates is 0(m). 
Notice that in recolouring a vertex v we might have changed the colours of neighbours 
of v as well. However, we never change the colour of a vertex whose colour agrees with 
the destination colouring a 2 , a fact that ensures that the process described above indeed 
terminates with the colouring a 2 . □ 

We are now able to show how to construct a flow such that A((f>) G 0(1) for q = 5 
colours. This only holds when the boundary colouring B of OR is the 0-colouring. 

Lemma 24. Suppose q = 5. Consider any region R and let B the the 0-colouring of OR. 
There is a flow such that the congestion A((j>) G 0(1). 
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Proof. For every pair (<7i,<7 2 ) G E d we know that o~\ and o" 2 differ only on vertices that 
are contained in the ball Ballot") for some vertex v G R d . Let -< be a fixed canonical 
ordering of the vertices in R. Let 7< Tl) < 72 G V ai , a2 be the path from o~\ to <x 2 constructed 
according to the proof of Lemma 23. We consider vertices in order specified by -< to 
make sure that 7< Tl)tJ - 2 is well defined. 

Assign all of the flow from U\ to a 2 to path 7< 71iCT2 G V ai ^ 2 . That is, 0(7 CT i,ct 2 ) = 
Prf(<Ji, cr 2 )7Tg((7i) and 0(7) = for all paths 7 G T^o^^o-i,^}- Let #1 and 9 2 , where 
(#1, 9 2 ) G £"0, be two colourings that disagree on a vertex w. Then the congestion of edge 
(0i A) is 

I T<T1 ,<T2 I ' " 



(<Tl,<T2)£E d : 
(Si,6 2 )G7o- 1 ,ct2 



WiA) 7r B (9i) 



^ E ».-^^ <*■•*.< 0(1), 

(o-l,cr 2 )eA d . 
(^l,^2)€7o- 1 ,(T2 

where fci and /c 2 are constants, specified next. Note that / kb{o~i)I'Kb{Oi) — 1- 

The path length |'y<Ti,<T 2 | * s upper-bounded by a constant k\ since <7i and <r 2 differ only 
on vertices inside a ball of fixed radius d. The path 7 CTli(T2 is constructed such that for 
each vertex v that is updated, we do at most a constant number of recolourings of vertices 
that are within constant distance from v. 

To see that the last sum is bounded by a constant k 2 , note that there are only a 
constant number of pairs (a"i,<7 2 ) in the summation. This is true since o~\ and <7 2 agree 
with 0\ on all vertices in R except in a constant-sized ball around a vertex w on which 
Q\ and 9 2 differ. Let m be the number of vertices u such that R d contains all vertices on 
which o~i and o 2 differ. Note that m is bounded by a constant since o~\ and a 2 differ only 
on vertices inside a ball of fixed radius d. We have 

1 



7TI i 

P^)< M eO( M ) 



Furthermore, 

Po((M2) >_L.I 6(i( _L ) 

since 1/g is the smallest probability of making a transition in Aio from colouring #1 to 
9 2 once vertex u> on which 9\ and # 2 differ has been chosen for an update. Thus, 

5M 6 o(i) 

Po(«iA) 1 J 

and we have that the sum is bounded by a constant k 2 . 

Now, Aq 1 o 2 ((f>) G 0(1) for all (6>i, 6> 2 ) G £"0 and it follows that the congestion A((f>) G 
0(1). " □ 
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Figure 9: The region Rx of the edge-boundary pair X in (a) is split into two halves (b) 
and (c). The split is along the vertices labelled v and vx- 

Finally we have the machinery for proving Theorem 5. 

Theorem (5, repeated). For any region R of the kagome lattice and q = 5 colours, the 
Glauber dynamics is rapidly mixing on R under the 0-colouring of OR. The mixing time 
t{8) G 0{n 2 + nlog |) ; where n is the number of vertices in R. 

Proof. The theorem is proved by using Lemmas 21 and 24. □ 

The proof of Theorem 6 is similar to the proof of Theorem 5. The implications from 
rapid mixing of to rapid mixing of the heat-bath dynamics on edges hold. Lemma 21 
has to be stated with A4o replaced by the heat-bath dynamics on edges (which slightly 
changes the proof) and Lemma 24 has be adjusted to deal with an arbitrary go-colouring 
of the boundary of the region, where q — 5. Showing that the congestion is constant 
under any g -colouring of the boundary is not difficult since we are allowed to update two 
vertices at the same time. 

8 The computational part of Lemma 15 

The computational part of the proof of Lemma 15 consists of two tasks: calculating the 
values f^F fm and constructing an (A, J 7 ) -set S that is good with respect to e — 1/1000. 
These two tasks are explained in the next sections. Both tasks are carried out using 
computer assistance. We have written programs in C, and the source code can be found 
on the webpage http : / /www . esc . liv . ac . uk/~markus/kagome5colours/ 

8.1 Computing jif^ m 

Calculating the values HF fm is a computationally challenging task. We are going to to 
calculate fJ>F fm for / £ {L • • • , 4720} and m G {1, . . . , 4}. From the definition of fJ-F fm in 
Section 5.3, ^F Um = if F' f m = 0. For every fixed / G {1, . . . , 4720}, F' f m = for exactly 
two values of m G {1, . . . , 4}. Thus, we will have to calculate the value of 2 x 4720 = 9440 
constants fJ>F fm - We must be able to compute a single value rather quickly, otherwise the 
total running time for all values will be too long. A brute-force approach would result in a 
running time of several months, maybe even years. We use a technique that is illustrated 
with the following example. 

Suppose TZp f is the extended region in Figure 9(a) and suppose m — 1. Hence the 
set F'f 7^ 0. The value (J>F fm is obtained by maximising fj,i j2 (X) over all edge-boundary 
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pairs X G F'f . Let v' be the vertex that is a neighbour to both vx and Since 
m — 1, the boundary edge between % and v' has colour 1 in every edge-boundary pair 
X G Ff m . An edge-boundary pair X G Fj m is therefore uniquely specified by the colour 
of its remaining boundary edges. Figure 9(a) illustrates an arbitrary edge-boundary pair 
X G F'f , where boundary edges are labelled with their colour (z\, . . . ,£13). Thus, in 
order to maximise H\ t2 {X) over all X G F'f , we could loop through all combinations of 
the colours zi, . . . , Z\% and compute fj,i i2 (X) for each such combination. This process will 
take very long. Next we explain how to speed up the process. 

By computing /i 1>2 (X) for many colourings of the boundary, one quickly makes the 
observation that only some particular colourings of the boundary result in a large value 
of /i lj2 (X). For other colourings, /i lj2 (X) tends to be rather small. For example, it turns 
out that if zi, . . . , zq are all colour 1, then /ii }2 (X) will be rather small regardless of the 
remaining colours 27, ... , z\%. Thus, setting the colours z±, . . . , zq to 1 is a "bad" choice if 
we want to maximise /ii t2 (X). From this observation we conclude that if we can filter out 
certain "bad" colourings of the boundary then we could speed up the process of finding 
the maximum value ^ 2 {X). 

We "split" the extended region 1Z Ff into two extended regions 1Z Mt and 1Z nght . Fig- 
ure 9(b) and (c) illustrate TZ Mt and TZ Tlght , respectively. The two regions share the vertices 
in the split. In this case it is vertex vx and v, both labelled in the figure. Let X lcft be 
the edge-boundary pair such that R x iat = 7Z Mt , e x ^ = ex and the boundary edges 
receive the same colours as in X. Boundary edges that are introduced from the split are 
given colour 0. Let X nght be the edge-boundary pair defined similarly to X lcft but with 
R Xlig u = n Tight . Figure 9(b) and (c) illustrate X left and X right . 

Let B be the colouring of £Rx such that B(e) = B x (e) for e G £Rx\{&x} and 
B(e x ) = 0. Recall from Section 5 that for i G Q, rii(X) denotes the number of proper 
g-colourings a in Qr x (B) such that u(vx) = i- For two colours i,i' G Q we now define 
n^°} h to be the number of proper g-colourings a in Q Rx (B) such that a(v x ) = i and 
a(v) = i', where v is the second vertex in the split. Thus, 



ni {x) = 



„both 
n i,i> ■ 



i'=l 



Let i? lcft be the colouring of £R X ^ such that B ieit (e x ^t) = and _B lcft (e) = B x ^t(e) 
for e G £R x ic{t\{e x ic{t}. For two colours i,i' G Q we define n- 6 ^ to be the number of 
proper g-colourings a in Q Rxlcft (B Mt ) such that cr(vx) = i and <r(v) = i' . We define n"f, ht 
similarly for the edge-boundary pair X right . It follows that 



^both _ left right 
u i,i> — u i,i' u i,i' 1 



and hence 



n t (X) = 



left right 
n i,i' n i,i' ■ 



= 1 
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With q = 5 colours, we have 



/ V N V^ 5 „ left „ right 



£.£{1,3,4,5} n * ( X ) Ei G {l,3,4,5} J2k=l n i?N"f ^ 



E 



„ left right 



y — — . (?) 



V Y^ nleft^rigti 

1 Z^ie{l,3,4,5} Z^fc=l n i,k n i,k 

1 



K , „lcft „ right 

J Z^e{i,3,4,5} 2^k=i U» x „right 

Note that the colours ^i, . . . , zq specify the quantity nfp , and z 7 , . . . , z± 3 specify the 

quantity n"f, ht . In order to maximise fii^(X) over edge-boundary pairs X, we could 
consider all combinations of the colours Zi,...,z±3 and use Equation (7). There are 
5 13 w 1.2 x 10 9 such combinations, so considering them all will take very long. Now, 
consider two different sets of the six colours Zi,...,z$. For G Q, let nfp' 1 be the 
value of nYP for the first set of colours, and let n) e p~ 2 be the value of n l fP for the second 
set of colours. Suppose 

n left-l n left-2 

„left-l — „ left- 2 \°) 
U l,j n l,3 

for all i G {1, 3, 4, 5}, j G {1, . . . , 5} and k G {1, . . . , 5}. Then we have from Equation (7) 
that /ii i2 (X) can only get smaller if we take nfp = n l fp~ 2 instead of nfp = nfp' 1 . In 
other words, there is no point considering the colours specified by the second set of 
colours zi,...,Zq when maximising /i 12 (X). This observation suggests that we loop 
through all combinations of colours Zi,...,z 6 and compare each pair of combinations 
like in Equation (8). We only keep the sets of colours that cannot be ruled out in some 
pairwise comparison like the second set above. This gives us a collection C lcft of colours 
Zi, . . . ,zq that turns out to be much smaller than the collection of all 5 6 = 15, 625 sets of 
colours. Similarly we obtain a collection C nght of colours z?, . . . , Zi 3 for the right part of 
the region. In order to find which colours z±, . . . , Z\ 3 that maximise /i li2 (X) we combine 
C left with C right . That is, we use Equation (7) to compute /ii^X) for each set zi, . . . , zq 
of colours in C lcft with each set 27, ... , z 13 of colours in C n s ht . 

The technique of splitting regions and filtering out boundary colourings that are guar- 
anteed not to maximise fii^iX) has a huge impact on the running time of the program. 
On a fairly powerful home-PC as of year 2006, it takes about two days to to obtain all 
9440 values a Ft ■ 



8.2 Constructing an (^4, jF)-set 

We describe how to construct an (^4, jF)-set. Let 7^big be the extended region in Fig- 
ure 10(a) with some combination of labels "in" and "out" on the vertices. From 1Z hlg we 
will derive 5-tuples that are added to a set S. By considering all possible combinations of 
labels "in" and "out" on the vertices of 7^ b i g) we construct the (^4, jF)-set S. We describe 
the process by first giving a concrete example. 
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Figure 10: (a) The extended region 7Z hig (here all vertices are labelled "in"). Note that 
the shaded vertex is not a part of the region, (b) The extended region TZa repeated. 
(c)-(f) Intersections of TZa and 7Zb ig . 

Fix an "in/out" -labelling of the vertices of the extended region TZ^g. Let a G 
{1, . . . ,342} be the value such that 7ZA a is an extended subregion of 7^bi g - Note that 
the extended regions TZa v i £ {1, • • ,342}, are defined such that there is exactly one 
value a G {1, ...,342} for which this is true. Figure 10(b) shows the largest possible 
lZA a and Figure 10(c) shows the overlapping of 7Zu g and 7Z,A a - We see from this figure 
that only some of the vertices of 7^big define 71 A a - Similarly to how the extended region 
7Za u is obtained from 7^big, let a±, a 2 , a 3 G {1, . . . , 342} be the three unique values such 
that r R-A ai is obtained from 7^big by the overlapping in Figure 10(d), 7ZA a2 is obtained 
from 7^-big by the overlapping in Figure 10(e), and 7Za uz is obtained from 7^big by the 
overlapping in Figure 10(f). It is possible that neighbours of vertex v in Figure 10(d)-(f) 
are labelled "out", meaning that some of the extended regions 7Z ai might not exist. If 
this is the case we define Oj = and A = 0. For example, if the vertex to the left of 
vertex v in Figure 10(d) is "out" then 7Z ai cannot exist and hence ai = 0. 

Suppose X is an edge-boundary pair such that Rx and 7^big are matching with respect 
toedgeex- Then X G A a . Fori G {1, 2, 3} and any two distinct colours G Q such that 
Px in (jjf) > 0, suppose Xi(j,f) is the extended edge-boundary pair that is constructed 
recursively in the tree Tx- If Xi(j,j') = then Xi(j,j') G A a .. We will now be more 
precise about the sets of edge-boundary pairs and incorporate the sets M 1; . . . , M 4 . 

Suppose without loss of generality that Bx{ex) = c and B' x (ex) = c', and /i CjC /(X) > 
/V jC (X) for some c, d G Q. Suppose the extended region in Figure 11(a) is an extended 
subregion of Suppose that the colour of the edge between wx and v in Figure 11(a) 
has colour c in Bx and B' x . Then X G Mi and hence X G A a ^. From Figure 11(a) 
we see that the extended region 7Zm (3 4) in Figure 4(b) is an extended subregion of 7ZA ai ■ 
Hence Xi(j,j') belongs to M 3 or M 4 (or both). The crucial observation here is that 
Px in (jif) > if and only if f = c. This follows from the fact that /i CjC /(X) > /i c / jC (X) 
and hence there is a discrepancy at vx only when the colour c is drawn from ttb' in the 
coupling \I/™ m . We therefore conclude that Xi(j,j f ) G M 3 . Thus, Xi(j,f) G A ai3 . For 
X?.{j,j') and X 3 (j,j') we see in Figure 11(a) that these edge-boundary pairs belong to 
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(3,0,0) (4,0,0) (0,0,0) 

Figure 11: Twelve cases which cover all possible combinations of the sets Mi, . . . , M4 to 
which the recursively constructed edge-boundary pairs Xi(ci, C2), X 2 (ci, c 2 ) and X 3 (ci, c 2 ) 
belong. If X E Mi then (a)-(d) apply. If X e M 2 then (e)-(h) apply. If X G M 3 U M 4 
then (i)-(l) apply. 
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either M 1 or M 2 . However, we are unable to tell exactly to which of the two sets these 
edge-boundary pairs belong. We therefore assume that any combination of the two sets 
is possible. The 3-tuples listed in Figure 11(a) indicate to which possible sets Mi, . . . , M 4 
the edge boundary pairs Xi(J,f), X 2 (j,f) and X 3 (j,j') belong. That is, a 3-tuple 
(m 1 ,m 2 ,m 3 ) means that Xi(j,f) £ A aumi , X 2 (j,f) £ A a2>m2 and X 3 (j,f) £ A a3>ms . 

Let T' C T be the set of edge-boundary pairs F i:1 such that 1Z Fi is an extended 
subregion of 7^big and i £ {1, . . . ,4720}. Then X £ F itl for every F it i £ J 7 '. Remember 
that we have assumed above that X £ Mi. Let / be the value such that Ff i £ T' is the 
set that minimises fj,p il over all F iti £ J 7 '. If the minimiser is not unique, let / be the 
smallest i among the minimisers. Now, for each 3-tuple (mi,m 2 ,m 3 ) in Figure 11(a) we 
add the following 5-tuple to the set S: (A aA , F fjl , A aumi , A a2jTn2 , A a2>m2 ). 

Summing it all up, we construct the set S as follows. First take an extended region 
IZbig- From Hug we uniquely derive the sets A a , A ai , A a2 and A az . If 1Zm {12) is an 
extended subregion of TZug then we consider two values of m: m — 1 and m — 2. If 
TZm, 3 4) is an extended subregion of Hug then we also consider two values of m: m = 3 
and m — 4. Now suppose X £ A a>m . The twelve cases in Figure 11 cover all possible 
combinations of the sets Mi, . . . , M 4 to which the recursively constructed edge-boundary 
pairs Xi(j,j f ), X 2 (j,j') and X 3 (j,j') belong. More precisely, if m = 1 then Figure 11(a)- 
(d) apply, if m = 2 then Figure ll(e)-(h) apply, if m = 3 or m = 4 then Figure 1 1 (i) — (1) 
apply. From TZug and the value of m, we uniquely derive the set Ff iTn to which X belongs. 
For each 3-tuple (mi,m 2 ,m 3 ) in the relevant case in Figure 11, we add the following 5- 
tuple to the set S: (A a , m , Ff jTn , A aumi ,A a2jTn2 , A a2)in2 ). If the value of m, in a 3-tuple is 
then A aumi = 0. By considering every possible extended region TZb lg and every possible 
value of m (two values per region TZhig), we construct a set S that is an (A, J 7 ) -set. 
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